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In a parameter regime for which the mean-field (Gross-Pitaevskii) dynamics becomes chaotic, 
mesoscopic quantum superpositions in phase space can occur in a double-well potential which is 
shaken periodically. For experimentally realistic initial states like the ground state of some 100 
atoms, the emergence of mesoscopic quantum superpositions in phase space is investigated numer- 
ically. It is shown to be reproducible even if the initial conditions slightly change. While the final 
state is not a perfect superposition of two distinct phase-states, the superposition is reached an 
order of magnitude faster than in the case of the collapse and revival phenomenon. Furthermore, a 
generator of entanglement is identified. 

PACS numbers: 03.75.Gg, 05.45.Mt, 03.75.Lm 

Keywords: double-well potential, mesoscopic entanglement, Bose-Einstein condensation 



I. INTRODUCTION 

Periodically shaking potentials used to trap ultra- 
cold atoms is currently established as an experimental 
method [IH3] : tunneling has been suppressed by choosing 
the correct amplitude of driving both for a Bose-Einstein 
condensate (BEC) in an optical lattice [H IH [5] and for 
a single atom in a double well potential O |6] . By fine- 
tuning the driving frequency, tunneling in tilted systems 
can also be enhanced in an analogue of photon-assisted 
tunneling. After the driving frequency has helped to 
bridge potential differences, tunneling can again be con- 
trolled by tuning the driving amplitude [U 2] [7] ■ Bose- 
Einstein condensates in periodically driven systems have 
been investigated regarding tunneling E] , dynamic lo- 
calization [5] and stabilities of exact Floquet states [10] . 

Suggestions how mesoscopic superpositions like the 
NOON-state QT] 

l*} = -L(|jv,o} + |o,iv}) (i) 

could be obtained with Bose-Einstein condensates can be 
found in Refs. [EIH25]; for a Bose-Einstein condensate in 
a double well this corresponds to a quantum superposi- 
tion of all N particles being either in one well or in the 
other. Here \n,N — n) denotes the Fock state with n par- 
ticles in the left well and N — n particles in the right well. 
Like the spin-squeezed states investigated in Refs. [27H29] 
such states are relevant to improve interferometric mea- 
surements [50] , 

A promising approach to obtain mesoscopic superpo- 
sitions are the collapse and revival phenomena in phase- 
space if the tunneling between both wells of a double- 
well potential is suppressed (or alternatively, for a single- 
species condensate in a harmonic trapping potential the 
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transition between two hyperfine-states) . A NOON-state 
then appears on the time scale [3T] 



where 2Hk is the interaction energy of a pair of atoms. 
Chaos-induced entanglement generation in a periodically 
shaken double well should, in principle, offer the possi- 
bility to obtain mesoscopic superpositions on short time 
scales [32] • While the state will not be a perfect super- 
position, the time scale on which it is reached will be an 
order of magnitude smaller. In this paper, the question 
is if these superpositions can be achieved both for ex- 
perimentally realistic conditions and for parameters such 
that the superposition is reproducible: in an experiment 
the particle number will vary from shot to shot and the 
tilt in a double well or the amplitude of shaking is likely 
to be slightly different. 

The paper is organized as follows: Section [Tl] intro- 
duces the model used to describe a Bose-Einstein con- 
densate in a periodically shaken double well. Section [TTI] 
describes experimental signatures of the entangled states. 
In Sec. |IV| a generator of mesoscopic superpositions is 
discussed while Sec. [V] demonstrates that the generation 
of mesoscopic quantum superpositions is robust against 
slight changes of the initial conditions. 



II. A BEC IN A PERIODICALLY SHAKEN 
DOUBLE WELL 

Bose-Einstein condensates in double-well potentials 
are interesting both experimentally and theoretically (T2J 
[28l 133- 37 . In order to describe a BEC in a double 
well, we use a model originally developed in nuclear 
physics |38| : a many-particle Hamiltonian in two-mode 



2 



approximation 39 , 



+ h(n + fix sin(wi)) (4c 2 - cjci) , 



(3) 



where the operator cy annihilates (creates) a boson in 
well j; Ml is the tunneling splitting, 2Hp,o is the tilt be- 
tween well 1 and well 2 and h^i is the driving amplitude. 
The interaction energy of a pair of particles in the same 
well is denoted by 2Hk. 

The Gross-Pitaevskii dynamics can be mapped to that 
of a nonrigid pendulum |34j . Including the term de- 
scribing the periodic shaking, the classical Hamiltonian 
is given by: 



H, 
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(4) 



where z is the population imbalance with z = 1 (z = — 1) 
referring to the situation with all particles in well 1 
(well 2). 

On the TV-particle quantum level, if all atoms occupy 
the single-particle state characterized by z = cos(6>) and 
<h. this leads to the wave-function 



JV 

\ 6 >ti = E(" ) cos"(0/2)sin w -"(0/2) 

n=0 ^ ' 

x e l(N - n)4, \n,N -n) . (5) 
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These bimodal phase-states [3T] are sometimes referred 
to as atomic coherent states. Note that for finite N these 
are not orthogonal: 



|(6» ,cj) \6\(j)')\ 2 > 0, N < oo. 



(6) 



As the Gross-Pitaevskii dynamics corresponds to that 
of a non-rigid pendulum, the dynamics is known to 
display a coexistence between chaotic and regular re- 
gions [40 . Here, the relation to entanglement [32] is 
discussed (cf. [JTJ 02] for the delta-kicked top) . 



III. SIGNATURES OF MESOSCOPIC 
QUANTUM SUPERPOSITIONS 

The signatures of entangled states described in this sec- 
tion can only serve as signatures of entangled states |43j if 
they are followed by a revival. However, as these revivals 
are only partial, this will not be a proof of entangle- 
ment but nevertheless a strong indication. Furthermore, 
on the level of computer simulations, one does not have 
the necessity to distinguish between quantum superpo- 
sitions and statistical mixtures. Deepening information 
concerning the states of interest, superpositions of atomic 
coherent states ([5|, can be found in [52"] . 



One signature of entanglement is the quantum Fisher 
information [44j [45] , here for the relative phase between 
the condensates in the two potential wells (in the fol- 
lowing abbreviatory referred to as the quantum Fisher 
information). Defining for pure states 



Fq F i = (An 12 ) 



(7) 



where (A7112) is the experimentally measurable [28] vari- 
ance of the particle number difference between both wells, 
one has an entanglement flag. For pure states, 



> 1, 



_ ^QFI 



N 



(8) 



is a sufficient condition for particle-entanglement and 
identifies those entangled states that are useful to over- 
come classical phase sensitivity in interferometric mea- 
surements [H] . Defining entangled states as states which 
cannot be written as product states [31] . this statement 
is supported by the fact that for a single atomic coherent 
state (5j) one has F ent < 1. This is expected since in this 
case all atoms occupy the same single-particle state. For 
mixtures, the quantum Fisher information becomes more 
complicated (Ref. [44] and references therein) and might 
not be easy to measure. Taking the definition (|7|) for an 
experiment leaves the difficulty to distinguish superposi- 
tions from mixtures. 

Another experimental signature [46] will be provided 
by interference patterns after switching off the potential 
and letting the wave-function expand for some time. A 
single atomic coherent state \9,4>) leads to an interference 
pattern used to experimentally detect the phase between 
the condensates [33] 



I(X) = [1 + sin(0) cos(A - 0)] 



(9) 



which still has to be multiplied by a Gaussian envelope; 
X essentially is the spatial variable in the direction con- 
necting the two wells of the double well and I the in- 
tensity at A. To characterize the quality of interference 
fringes, the visibility or contrast 



C = 



+ In 



(10) 



is introduced where / max is the highest and / m i n the low- 
est value of /. Given the fact that even a single atomic co- 
herent state can produce an impressive interference pat- 
tern, interference cannot be used in a straightforward 
manner as an entanglement flag as in Refs. [23J[3B]. Let 
us thus consider the disappearance of interference pat- 
terns, 



C<1 



(11) 



This might simply mean sin(0) ~ (cf. Eq. Q) or even 
heating. However, when combined with a high quantum 
Fisher information, it will be the signature for the super- 
position of two (or more) distinct atomic coherent states 
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FIG. 1: (Color online) Comparison of Poincare surface of 
section (a) and QFI-map (b) for parameters 2/j,o/£l = 3.0, 
Nk/VI = 0.5, w/f2 = 3.0 and 2/ii/w = 0.1. Plot (b) shows the 
time averaged QFI for a system with N = 100 particles and 
averaging time 10/O. The QFI is enhanced close to the hyper- 
bolic fixed point in the Poincare section of surface (z w 0.21, 
cj> « 0.54-7T), which acts as a generator of mesoscopic entan- 
glement (c). 



which are only present in mean-field chaos, at least some 
of the aspects should be understandable in a parameter 
regime for which the mean-field dynamics is regular and 
thus the corresponding Poincare surface of section (cf. 
Fig. [l^) is easy to understand: essentially, there is one 
elliptic fixed point at z m 0.17, <j> » — 0.487T and one 
hyperbolic fixed point near z 0.21 and <p ~ 0.547T, 
in addition there are further fixed points near the top 
and bottom of the plot which can be discarded at this 
point. Within the mean-field dynamics, a solution which 
starts right on a hyperbolic fixed point (#h yP cr, 0hyper) 
will not move at all. This, however, is different on the 
TV-particle quantum level: preparing the system in the 
atomic coherent state corresponding to all N particles 
being in the atomic coherent state |#hypen </>hyper) will - 
for any N < oo - have contributions from neighboring 
atomic coherent states, cf. Eq. ([6]). Those states, even 
on the mean-field level, do move. 

Figure [TJd shows the time-averaged quantum Fisher in- 
formation of all particles initially occupying an atomic 
coherent state. For starting values in the vicinity of the 
hyperbolic fixed point mentioned above, mesoscopic su- 
perpositions occur. If this is indeed what is happening, 
the remaining question is why the other hyperbolic fixed 
points do not display such a behavior: contrary to the 
first one, they do not lie near the separatrix for which 
the dynamics happens on much shorter time scales. 

Related to the observed acceleration of entanglement 
generation is the identification of enhanced atomic tun- 
neling between the two wells in high-chaoticity regions 



with sin(#i) ~ sin(# 2 ) and cos(X — (f>i)+cos(X — <j) 2 ) = 0. 
This is a strong indication for mesoscopic superpositions 
similar to the case of the superfluid to Mott-insulator 
transition [37] where the dis- and reappearance of the 
contrast clearly indicates (but does not proof) the Mott- 
insulator transition. 

In order to numerically calculate the contrast, one 
needs 1351 



with 



\¥(x)^(x)\i;) 



*(.t) = $i(a;)ci + $ 2 (a;)c2, 



(12) 



(13) 



where "3^ is the expanded mode which was originally lo- 
calized in well i with i = 1,2 before switching off the 
potential. 



IV. EMERGENCE OF QUANTUM 
SUPERPOSITIONS 

The claim of Rcf. [32] is not that mean-field chaos is 
the origin of entanglement generation but rather that it 
speeds up the emergence of mesoscopic superpositions. 
Thus, rather than looking for entanglement generators 



V. ROBUSTNESS OF ENTANGLEMENT 
GENERATION 

In this section parameters are presented where an en- 
tangled state occurs on short time scales, indicated by 
both signatures quantum Fisher information and con- 
trast. The robustness of this state under realistic experi- 
mental conditions is discussed. In the following the initial 
state is prepared without driving, hence corresponds to 
the thermal occupation of the eigenstates with energies 
Ei of Hamiltonian (3| exclusive of the last term. In a 
typical experimental situation only some energy eigen- 
states are significantly occupied 28 which could consid- 
erably be reduced by working at lower temperatures in 
the sub nano Kelvin range [SO]. The initial state can be 
prepared in a well-controlled manner and the experimen- 
tal parameters are well controllable [55J except for small 
experimental uncertainties. The occurring fluctuations 
of the initial particle number are discussed in this sec- 
tion while particle losses are not analyzed; their effect on 
phase revivals has been investigated in [51) . 

Figure [2] shows the time development of quantum 
Fisher information and contrast for initially ground and 
first excited state in dependence of the scaled driving am- 
plitude 2fii/u. Obviously both initial states show a sim- 
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FIG. 2: (Color online) Entanglement flag F cnt as a func- 
tion of time t — Qt and scaled driving amplitude 2/xi/cj for 
ground (a) and first excited state (b). (c) and (d): Contrast 
for ground and first excited state. N = 110, Nk/Q = 1.045, 
fio/Q, = 0.75, w/fi = 3.0. 




FIG. 3: Entanglement flag -F cnt (a) and contrast (b) for 
2/ii/o; = 1.0. Other parameters as in Fig. [5] 



ilar behavior. Parameters where both quantum Fisher 
information and contrast strongly indicate entanglement 
are shown in Fig. [3] The quantum Fisher information 
takes on a value of 70.5 and the contrast decreases to 
1 % already at r = 7.1. Additionally the contrast reap- 
pears a short time later, hence heating can be excluded 
in an experiment. 

To determine quantum Fisher information or contrast 
in an experiment it has to be repeated several times with 
the same initial conditions. But in the laboratory param- 
eters like driving frequency uj and tilt 2hfj, will change 
slightly from shot to shot. Figure [4] shows the influence 
of slight deviations on the two entanglement indicators. 
Another parameter where deviations will occur is the ini- 
tial particle number N. The dependence of the signatures 
is represented in Figs. [5] For Poissonian distributed ini- 
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FIG. 4: (Color online) Entanglement flag (a) and contrast (b) 
at r = 7.1 versus fi /fl and lu/Q. N = 110, Nk/Q = 1.045, 
2/iti/w = 1.0. 



tial particle numbers the signatures are likely to be still 
visible in the experiment. 

So far, we have optimized with respect to interfer- 
ence patterns and fluctuations of the initial particle num- 
ber N. Other experimentally interesting quantum states 
are characterized by a bimodal phase-distribution with a 
small standard deviation. Using the completeness rela- 
tion [52] 



N+ 1 

47T 



d0sin(0) / d0| 



0| (14) 



the phase distribution of a state can be calculated. 
The initial particle number is again assumed to obey a 
Poissonian distribution with mean value N = 110 which 
for numerical purposes has been truncated and renormal- 
ized with N ranging from TV = 100 to N = 120. Then, a 
mean phase distribution and its standard deviation can 
be calculated for each point of time which has been opti- 
mized in order to find a phase distribution that indicates 
strong entanglement and as well has a small standard de- 
viation. In such an optimization procedure an entangled 
state has been found at r = 6.5 where the phase distribu- 
tion shows two distinct maxima, separated by A<f) ir, 
cf. Fig. [6]d. For maximum reappearance one important 
peak is observed at r = 8.9 (c). For comparison, an ex- 
emplary mixed state leading to the phase distribution of 
(b) has been constructed: the wavefunction at r = 6.5 
has been split into two parts, each consisting of the con- 
tributions to Prj, in an interval of A(f> = tt around one of 
the peaks. In the statistical mixture these two, renormal- 
ized parts are weighted according to their contribution to 
the total wavefunction. This mixed state would after the 
corresponding time difference At = 2.4 have a phase dis- 
tribution according to (d) which is distinct of (c). Hence, 
it is possible to distinguish pure entangled states from 
mixed states. Additionally, preliminary results indicate 
that under slight parameter variations of about 1 % for 
driving frequency, driving strength and initial tilt respec- 
tively 10 % for the scaled interaction strength Nk/SI, the 
phase distribution at r = 6.5 and the partial revival at 
r = 8.9 stay nonetheless robust. 

Besides slight changes in the experimental conditions 
another effect that should be discussed is the influence of 
temperature. So far, calculations have been performed 
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FIG. 5: (Color online) Font (a) and contrast (b) for r = 7.1 
(black crosses) and r = 9.0 (blue stars) for different particle 
numbers such that /t/Q = const. Nk/Q(N = 110) = 1.045, 
Ho/0, = 0.75, w/fi = 3.0, 2^i = 1.0. 
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FIG. 6: (Color online) Averaged phase distribution p$ 
(red solid line) with standard deviation (blue dashed lines). 
The initial particle number is modeled by a truncated, 
renormalized Poisson distribution of the initial particle 
number, ranging from N = 100 to N = 120, with 
mean value N = 110. (a) Initial state, (b) entangled state 
at t = 6.5, (c) time-developed state at r = 8.9, (d) 
assumed exemplary, two-component statistical mixture at 
r = 8.9. The interaction strength n/Q. is kept constant 
such that Nk/Q(N = 110) = 1.0. fi /n = 1.3, u/Q. = 3.3, 
2/n/oj = 0.7. 



for T = 0. To account for temperature the eigenstates 
of the undriven system corresponding to energies Ei are 
supposed to be populated with probability pi according 
to a Boltzmann distribution 

p -(Ei-E )/k B T m 

Pi = j , Z = J2 e-^- E °y** T . (15) 

Hence, the statistically averaged contrast is given by 



M 
z=0 



(16) 



where Cj is the contrast for the pure state i and M is the 
number of states included in the average. With an ana- 
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FIG. 7: (Color online) F cnt for r = 7.1 (a) and r = 9.0 (c). 
Black, solid line: quantum Fisher information -Fent.mixed f° r 
mixed states (cf. Eq. (18l), blue, dashed line: {F cnt )st statis- 
tically averaged over all 111 states, red, dotted line: (-F e nt)st 
statistically averaged over three lowest lying, experimentally 
important states [28] ■ (b ) and (d): contrast (C) s t for r = 7.1 
and t — 9.0 (cf. Eq. (16l). Blue, dash-dotted line: statistical 
average over all states, red, dotted line: statistical average 
over three lowest lying eigenstates. Parameters same as in 
Fig. 3, 2/ir/w = 1.0. 



logue definition of (-) st the statistically averaged quan- 
tum Fisher information becomes 

/J? \ ((A n 12) 2 )st 1 ,. 2 x / \2\ / 17 \ 

(-Pent)st = ^ = ^ (K 2 )st - MrtJ ■ ( 17 ) 

For comparison also the quantum Fisher information for 
mixed states (Ref [H] and references therein) 



cnt, mixed 



2N ^ Pj + Pk 

3i K 



again relative to the relative phase, is calculated. Tem- 
perature effects are represented in Fig. [7] for the same 
parameters as in Fig. 2 and the interesting points of 
time r = 7.1 and r = 9.0 where the contrast vanishes 
respectively reappears maximally. As the experiments in 
|28j are performed in a temperature regime where only 
the three lowest lying states are significantly populated 
both the statistical averages over the three lowest lying 
states and all states are calculated. The curves for only 
the three lowest lying and all states stay almost identical 
up to temperatures where k^T approximately equals the 
energy separation of ground and first excited state. The 
value of the quantum Fisher information at r = 7.1 is still 
adequate while the contrast stays low. Also for r = 9.0 
the signatures are not influenced too much. For higher 
temperatures the quantum Fisher information for mixed 
states deviates from the statistically averaged quantum 
Fisher information (17). Hence, it should be noted that 



measuring the variance of the particle number difference 
at higher temperatures is not equivalent to measuring 
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the quantum Fisher information for mixed states while 
for low enough temperatures both quantities lead to sim- 
ilar results. 



Conclusion 

To summarize, the occurrence of entanglement on 
short time scales is observed for a Bose-Einstein conden- 
sate in a tilted double-well potential. For this system on 
the mean-field level the entanglement generation is ac- 
celerated in the vicinity of hyperbolic fixed points. The 
robustness of the entangled states under realistic exper- 
imental conditions is analyzed and parameters are pre- 
sented where the signatures quantum Fisher information 
for the relative phase and contrast are found to be suffi- 
ciently stable. In conclusion, for the proposed parameters 
the phase distribution of the entangled state measured in 



an experiment is expected to exhibit two clearly distin- 
guishable maxima. 

Our results are based upon the widely-used two-mode 
approximation. An analysis of its limitations can be 
found, e.g., in [S3]. Nevertheless, the two- mode approx- 
imation has, especially for comparatively small interac- 
tion strengths, been found to describe many aspects of 
the experiments [28, 33]. Besides the investigation of ex- 
tended models, the influence of particle losses is certainly 
an interesting aspect for future work. 
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